Counter-ions at Charged Walls: Two Dimensional Systems 
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We study equilibrium statistical mechanics of classical point counter-ions, formulated on 2D Eu- 
clidean space with logarithmic Coulomb interactions (infinite number of particles) or on the cylinder 
surface (finite particle numbers), in the vicinity of a single uniformly charged line (one single double- 
layer), or between two such lines (interacting double- layers). The weak-coupling Poisson-Boltzmann 
theory, which applies when the coupling constant F is small, is briefly recapitulated (the coupling 
constant is defined as F = fie 2 where /3 is the inverse temperature, and e the counter- ion charge). 
The opposite strong-coupling limit (r — !> oo) is treated by using a recent method based on an exact 
expansion around the ground-state Wigner crystal of counter-ions. The weak- and strong-coupling 
theories are compared at intermediary values of the coupling constant V — 2"/ (-y — 1, 2, 3), to exact 
results derived within a ID lattice representation of 2D Coulomb systems in terms of anti-commuting 
field variables. The models (density profile, pressure) are solved exactly for any particles numbers 
iV at r = 2 and up to relatively large finite N at F = 4 and 6. For the one- line geometry, the 
decay of the density profile at asymptotic distance from the line undergoes a fundamental change 
with respect to the mean-field behavior at V = 6. The like-charge attraction regime, possible in the 
strong coupling limit but precluded at mean-field level, survives for r = 4 and 6, but disappears at 

r = 2. 
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PACS numbers: 82.70.-y, 82.45.-h,61.20.< 



I. INTRODUCTION 

Most mesoscopic objects, when dissolved in a polar 
solvent such as water, acquire an electric charge through 
the dissociation of functional surface groups [l[ . Counter- 
ions are then released in the solution, and form, together 
with the charged object, the so-called electric double- 
layer. Since the pioneering work of Gouy and Chapman 
a century ago Q , the study of these charge density clouds 
has formed an active line of research, in particular from 
a theoretical perspective . Electric double layers are 
indeed pivotal in affecting single mesoscopic "particle" 
properties, together with inter-particle interactions. 

The present paper concerns the equilibrium statistical 
mechanics of charged particles in the vicinity of charged 
walls (planar double-layers). The general problem of mo- 
bile ions confined by uniformly charged interfaces can be 
formulated in two ways. In the case of "counter-ions 
only", there is just one species of equally charged ions 
neutralizing the surface charge on walls. In the case of 
"added electrolyte" , the system is in contact with an in- 
finite reservoir of ± electrolyte charges. Here, we shall 
restrict ourselves to the former class of models (no salt 
case). The mostly studied 3D geometries are one planar 
wall with counter-ions localized in the complementary 
half-space and two parallel planar walls with counter- 
ions being localized in between. The high-temperature 
(weak-coupling) limit of general Coulomb systems is de- 
scribed by the Poisson-Boltzmann (PB) mean-field the- 
ory [?} which basically does not depend on dimension. 
Formulating the Coulomb system as a field theory, the 
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PB equation can be viewed as a first-order term of a sys- 
tematic expansion in loops [sl-[TT|. 

A relevant progress has been made in the last decade 
in the opposite low-temperature (strong-coupling, SC) 
limit |p, ll2H2a |. Within a field-theoretical treatment 
[3 Il6j . the leading SC behavior stems from a single- 
particle picture and next correction orders correspond to 
a virial/fugacity expansion in inverse powers of the cou- 
pling constant. The method requires a renormalization 
of infrared divergences via the electroneutrality condi- 
tion. A comparison with the Monte-Carlo simulations 
fl7j confirmed the adequacy of the leading single-particle 
theory, but the predictions for the first correction turns 
out to be incorrect. In essence, a virial-like expansion 
fails, as it also fails for simple electrolytes (where a rem- 
iniscent divergence of the second virial coefficient is in- 
dicative of a non analyticity in the density expansion of 
the pressure, and signals that an alternative theoretical 
route should be explored [HI)- Very recently, a method 
based on an exact expansion around the ground state 
(the 2D Wigner crystal of counter-ions formed on the 
surface of charged walls) was proposed in Ref. [13] to 
overcome such shortcoming. The expansion is systematic 
and free of divergences, it provides the correct correction 
to leading SC behaviour, and the obtained results are in 
excellent agreement with available data of Monte-Carlo 
simulations under strong couplings [27| . 

The majority of previous studies has been restricted to 
3D which is the dimension of practical interest. Among 
notable exceptions are references |28|,[29] , dealing with 2D 
Coulomb systems of counter-ions with logarithmic pair- 
wise interactions, mainly in the context of Manning con- 
densation. Here, we shall concentrate on 2D systems of 
counter-ions, with charged lines on boundary walls, and 
Ref. [23d in one dimension. The homogeneous line charge 
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dition Vw(r) —> as r —> oo, reads 





FIG. 1. Considered geometries: a) one infinite charged line; 
b) two parallel charged lines; c) the semi-infinite cylinder of 
circumference W; d) finite cylinder. Hatched domains are 
accessible to particles, ae is the uniform charge density of 
lines. 



density is ere (e is the elementary charge and a > 0); 
point-like counter-ions are for simplicity monovalent with 
charge — e. The studied geometries are pictured in Fig. 
1; white walls are impenetrable to particles, domains ac- 
cessible to particles are hatched. The geometries of one 
(infinite) charged line and two parallel charged lines at 
distance d are presented in Figs, la and lb, respectively. 
The semi-infinite and finite cylinder surface area of cir- 
cumference W with charged circle boundaries are pic- 
tured in Figs, lc and Id, respectively. These cylindric 
geometries will enable us to mimic the former ones, ob- 
tained as the infinite-particle-number limit Wa — > oo, for 
finite numbers of counter-ions. The walls and the confin- 
ing domains are assumed to possess, for simplicity, the 
same (vacuum) dielectric constant e — 1, so that there are 
no image charges. The relevant dimensionless coupling 
constant is T = /3e 2 , where ft is the inverse temperature. 

Interactions of counter-ions with each other and with 
charged surfaces are determined by 2D electrostatics de- 
fined as follows. In v spatial dimensions, the electrostatic 
potential v at a point r € R", induced by a unit charge 
at the origin 0, is the solution of Poisson's equation 



Av(r) 



A*), 



(i.i) 



where s v is the surface area of the ^-dimensional unit 
sphere; S2 — 27r, S3 = 47r, etc. This definition of 
the ^-dimensional Coulomb potential maintains generic 
properties -such as screening sum rules- of "real" 3D 
Coulomb systems with the interaction potential v(r) — 
1/r, r g R 3 and r = |r|. In an infinite 2D Euclidean 
space, the solution of (jl.lj) . subject to the boundary con- 



v(r) 



In 



r G R 2 , 



(1.2) 



where the free-length scale L will be set for simplicity to 
unity. Such potential is created by infinitely long charged 
lines in 3D which are perpendicular to the given 2D plane; 
the corresponding systems occurring in the nature are 
the so-called polyelectrolytes. The Coulomb potential 
(|1.2[) is used for the geometries in Figs, la and lb. For 
cylindric geometries in Figs, lc and Id, the requirement 
of periodicity along the y-axis, with period W, leads to 
the Coulomb potential 30 



v(r) = — In 



2 sinh — 



(1.3) 



where the complex notation z = x + iy is used for point 
r = (x,y). At small distances r <C W, this potential 
behaves like the logarithmic one (ll.2p with L = W/(2tt). 
At large distances along the cylinder x ^> W, it behaves 
like the ID Coulomb potential —{tt/W)\x\. 

Maintaining basic features of physical phenomena, 2D 
models with one type of charges have interesting advan- 
tages in comparison with 3D ones: They are less labori- 
ous and some of the concepts can be often verified by ex- 
plicit calculations. 2D Coulomb systems are even exactly 
solvable at a special coupling T = 2, in infinite space as 
well as in inhomogeneous semi-infinite or finite domains; 
for a review, see [3l|. Such exactly solvable models can 
serve as an adequacy test of weak- and strong-coupling 
series expansions for finite temperatures. It was shown 
in Ref. 32| that for the sequence of couplings T = 2j 
(7 = 1, 2, 3, . . .) statistical averages in 2D Coulomb mod- 
els can be treated within a ID lattice theory of inter- 
acting fields of anti-commuting (Grassmann) variables. 
This ID representation enables one to treat exactly one- 
component Coulomb systems with relatively large num- 
bers of particles also for T — 4, 6. The technique of anti- 
commuting variables was applied to the general problem 
of integrability of the 2D jellium [33| and to the trans- 
lational symmetry breaking of the jellium formulated on 
the cylinder surface [34|. 

The present work was motivated by the need to have 
a control over weak- and strong-coupling theories. The 
exact results for density profiles and pressures at inter- 
mediary values of F provide valuable tests of weak- and 
strong-coupling theories. Among the results obtained in 
this paper, two deserve special attention. For the one- 
line geometry, the asymptotic decay of the density pro- 
file from the line undergoes a fundamental change from 
the mean- field behavior at T = 6. This means that 
the long-distance predictions of the PB theory have a 
restricted validity, which is relevant from the point of 
view of the renormalized-charge concept (for a review, 
see [Hj])- This also invalidates the hypothesis that a 
strongly coupled double-layer, at large distances, behaves 
in a (suitably renormalized) mean- field fashion [35l [36j j . 
For two-line geometry, there is evidence about attraction 



3 



between like-charged lines at relatively small couplings 
r = 4, 6. The attraction is observed even in the N = 2 
particle systems. 

The paper is organized as follows. For completeness, 
it starts with a short recapitulation of the weak-coupling 
PB theory for the straight-line(s) geometries in Figs, la 
and lb. The models are subsequently studied in the SC 
limit, by using the new method [27], in Sec. 3. Inter- 
mediary values of the coupling constant T are investi- 
gated within the technique of ID Grassmann variables in 
Sec. 4. The models, formulated on the cylinder surface 
in Figs, lc and Id, are solved exactly (density profile, 
pressure) for any particles numbers N at T = 2 and for 
(relatively large) finite N at T = 4, 6. The results are 
compared with those obtained in the weak- and strong- 
coupling limits. Conclusion are finally drawn in section 
5. 



II. WEAK-COUPLING LIMIT 

We begin with a brief reminder of the weak-coupling 
PB theory, adapted for the geometries pictured in Figs, 
la and lb. More details can be found in e.g. Q- 



A. Single charged line 

We first consider the case of a single infinite line at 
x = 0, carrying positive charge density ere (Fig. la). 
Let the density of counter-ions in the half-plane x > 
be denoted by n(x); the corresponding charge density is 
p{x) = —en(x). The contact theorem for planar wall 
surfaces |37H39j | relates the total contact density of par- 
ticles to the surface charge density on the wall and the 
bulk pressure of the fluid P. For 2D systems of identical 
particles, it reads 



PP = n(0) - irTa 2 



(2.1) 



where T = (3e 2 is the coupling constant. Since in the 
present case of a single isolated double-layer, the pressure 
vanishes (P — 0), we have n(0) = irTa 2 . 

The induced average electrostatic potential (f>(x) is de- 
termined by the Poisson equation 



d 2 <t>(x) 
da; 2 



= —2irp(x) 



(2.2) 



The condition of electroneutrality ere + dx p(a 
is equivalent to the boundary condition 



d(j)(x) 



dx 



2nae. 



x=0 



(2.3) 



Within the mean- field approach, valid for T — > ("high 
temperatures"), the average particle density is approxi- 
mated by replacing the potential of mean force by the 
average electrostatic potential, n(x) = uq exp[/3e0(a;)]. 



Equation (|2.2[) then reduces to the nonlinear PB equa- 
tion 



d 2 <j)(x) 
da; 2 



= 2neno exp[/3e<fi(x)} 



(2.4) 



This second-order differential equation, supplemented by 
the boundary condition (|2.3p . can be integrated explic- 
itly. The solution for the counter-ion density profile reads 



n(x) 



1 



1 



7rr (a; + /i) 2 



(2.5) 



Here, fi = l/(itTa) is the Gouy-Chapman length, i.e. 
the distance from the charged wall at which an isolated 
countcr-ion has potential energy equal to thermal energy 
k^T = 1//3. In what follows, all lengths will be expressed 
in units of fi, x — x/p. Note that at asymptotically large 
distances from the wall x — > oo, the density profile (|2.5|) 
does not depend on the line charge density magnitude 
ere, n(x) ~ l/(7rra; 2 ). We shall see that this interesting 
phenomenon takes place also at finite temperature (r = 2 
and r = 4), while another asymptotic decay seems to 
hold at T = 6. 

The density profile will be considered also in the 
rescaled form h(x) = n(fj,x) /(ttT a 2 ). Thus, the solution 
(|2.5[) can be expressed as 



h{x) = 



1 



1 



(1 



(2.6) 



This density profile satisfies both the electroneutrality 
condition (|2.3p . written in the rescaled form as 



dxn(x) = 1, 



(2.7) 



and the contact theorem (12.11) with P = 0, written as 



n(0) = 1. 



(2.8) 



It is worth emphasising that as an approximate approach, 
PB could a priori violate the contact theorem. This is 
not the follows from arguments that can be found 

in Ref. O. 



B. Two charged lines 

The system of two lines at distance d in Fig. lb has 
the x — ¥ d — x symmetry, so it is sufficient to consider 
the interval x e [0, d/2}. Because of the mentioned sym- 
metry, the electrostatic potential satisfies the additional 
boundary condition at the mid-point x = d/2 



d(j){x) 



dx 



x=d/2 



0. 



(2.9) 



Let us denote by n m the value of the particle density 
at x — d/2. The PB equation (|2.2I) . supplemented with 
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the boundary condition (|2.9[) . can be solved explicitly: 



/3e<p(x) = — In cos 2 
n(x) — 



K\x- 



[K(x - d/2)] ' 



(2.10) 



where the inverse length K is related to n m via 

K 2 = nTn m . (2.11) 



The boundary condition at x = (|2.3p implies the fol- 
lowing transcendental equation for the real K 



{Kd) t&u{Kd/2) = - = d. 



(2.12) 



According to the contact theorem (|2.1I) . the (rescaled) 
pressure between the charged lines is given by 



P = 



J3P_ 
nTa 2 



(2.13) 



Note that in the PB approximation, the pressure between 
two equivalently charged lines is always positive. This 
repulsive behaviour is in agreement with general results 
41]. 

Two limiting cases of the dimensionless distance be- 
tween the lines d are of special interest. For small dis- 
tances d — > 0, we have 



d 3 



while for large distances d 

f3Pr 



Id- 

45 



0(d 2 ), 



> oo, Kd —> 7r and 

7T 1 

fd 2 ' 



(2.14) 



(2.15) 



We see that the asymptotic decay of P to zero does not 
depend on a; this behavior is analogous to that of the 
density profile in the one-line problem. Note that d and 
l/cr (or equivalently, the Gouy length) being the two rel- 
evant length scales in the problem, an algebraic decay of 
P in 1 /d 2 at large distances cannot involve any a depen- 
dence, for dimensional reasons. A similar remark holds 
for the density profile of the single plate problem, Eq. 
(USD. 



a is much larger than the characteristic Gouy distance 
[x between the counter- ions and the charged line [Hj], 
a = a/fx oc r — > oo. In the asymptotic SC limit T — > oo, 
each vertex Rj is occupied by a counter- ion j. The 
ground-state energy of the discrete particle system plus 
the homogeneous background charge density will be de- 
noted by Eq. For T large but finite, the fluctuations of 
counter-ions around their vertex positions become im- 
portant. Our approach consists in a systematic account 
of these fluctuations, and provides an expansion in in- 
verse powers of r. It is the exact adaptation to the 2D 
case of the procedure discussed in 27| for three dimen- 
sional systems. It differs from the approach of Netz and 
collaborators [f| [IH, [I?} in that the latter virial-like pro- 
cedure, although capturing the leading order term in the 
expansion, leads to divergences in the corrections. As 
a consequence, the corrections to the leading order be- 
haviour are not of the order in the coupling constant that 
is predicted in [|| [T(| [l7j . see [27| ■ On the other hand, 
our expansion, that purports to capture the very same 
phenomena as Refs [g, [lfs, [l7|, is free of divergences, and 
reveals that the corrections to the leading order are much 
larger than obtained in [f| [l6|, [l7|- This is confirmed 
by the Monte Carlo simulations, that are, for three di- 
mensional systems, in complete agreement with the pre- 
dictions of [271 ] . Our strong coupling expansion method 
also differs from that put forward in refs [pj, [PJ] by Lau, 
Pincus, and collaborators. These authors impose that 
the counter-ions stick to the charged interfaces, so that 
the counter-ions degrees of freedom are "in-plane" only 
(along the interface). This precludes the possibility to 
study the phenomenon of like-charge attraction at small 
distances, see below. Indeed, it is essential to include in 
the analysis the excitations where the counter-ions can 
unbind from the charged interface (displacements per- 
pendicular to the interface). Consequently, the results of 
Refs. [TH, [l9| should be viewed as a large distance ex- 
pansion, and in this respect, complementary to our short 
distance analysis. 

As depicted in Fig. 2. a, let us first shift one of the 
particles, say j = 0, from its lattice position Ro = (0, 0) 
by a small vector 5R — (x > 0,y), SR <C a- The 
corresponding change in the total energy 6E(x,y) = 
E(x, y) — Eq > consists of two contributions. The first 
one is due to the interaction of the shifted counter-ion 
with the uniform line charge density: 



III. STRONG-COUPLING THEORY 

A. Single charged line 

At zero temperature, i.e. in the strong coupling (SC) 
limit r — > oo, the counter-ions collapse on the charged 
line, see Fig. 2a. They create a ID Wigner crystal, 
with vertices R 3 = (0,j/a) (j = 0, ±1, ±2, . . .); the 
nearest- neighbor distance a = l/cr ensures the local neu- 
trality on the line. It is essential here to bear in mind 
that the SC limit corresponds to the regime in which 



(5£ (1) (a 



7re ox. 



(3.1) 



The second contribution SE^ 2 '(x,y) comes from the in- 
teraction of the shifted particle with all other ions on the 
ID Wigner crystal. For the special case of the y = 
shift, we find 



5E^(x,0) = -e 2 J2 




In 
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-e <> 



-e 



a) 



1/a 



b) where use was made of ^2jLi 1/j 2 = t>" 2 /6 [42]. Note the 
absence of the mixed harmonic term xy. The total energy 
change, up to harmonic terms and in dimensionless form, 
is finally given by 



FIG. 2. Ground-state Wigner crystal of counter-ions (shown 
by the black dots) in a) one-line geometry; b) two-lines ge- 
ometry. Figure a) also introduces the notation (x,y) for the 
position shift of a tagged counter-ion, that will be used in our 
analytical calculation. 



mn 

3=1 



ox 

J 



s 2 In 



sinh(7rcra;) 



(3.2) 



This function has a small- x expansion of the form 

2 2 

5EW(x,0) = - — (ttctx) 2 + —(Trax) 4 + 0(x 6 ). (3.3) 
6 180 

Note that the interaction with counter-ions does not im- 
ply any contribution linear in x. The corresponding elec- 
tric field indeed vanishes, by symmetry. The negative 
sign of the leading x 2 term does not represent any prob- 
lem: The sum of 8E^ (x, 0) and the linear term (|3.1[) is 
a monotonously increasing function of x, as it should be. 
For the special case of the x = shift, we find 



5E^(0,y) = -e 2 lnJ] 



1-1^ 



-e 2 ln 



3=1 

sin(iray) 



Tray 



— {■Kayf + YgQ^y) 



C% 6 ).(3.4) 



The calculation of the whole function 5E^ {x, y) for the 
particle shift simultaneously along both directions is com- 
plicated. For our purpose, it is sufficient to derive its 
expansion in x,y up to harmonic terms: 



e 2 ^ 



-i 2 



6E^(x,y) = -yl> 

3 = 1 



<yy_ 
3 




l 



PSE(x,y) x + — (i 

61 



(3.6) 



This formula reveals a relationship between the order of 
the expansion of —/3E(x, y) in the dimensionless lengths 
x,y and its SC expansion in 1/r. The linear term —x is 
the only one which does not vanish in the limit T — > oo. 
This leading term reflects the single-particle character of 
the system close to zero temperature: Each particle be- 
haves independently from the others, exposed only to the 
field of the line charge density. A similar conclusion is 
reached within the original approach of Netz and collab- 
orators @, [l6l Terms of the pth order in x, y have 
prefactors proportional to Ta p , i.e. passing to x, y they 
become of the SC order l/r^" 1 . This is why the first 1/r 
correction to the single-particle regime is determined ex- 
clusively by the terms harmonic in space. 

Let us now consider a shift of all particles from their 
lattice positions Rj (j — 0, ±1, ±2, . . .) by small vectors 
<5Rj = (Xj , yj ) . To determine the corresponding increase 
in the total energy 5E({xj, yj}), we proceed as above and 
obtain 



p5E({ Xj , yj })~-J2- 



1 



27r 2 r 



i 



27r 2 r 



E 

j<k 



E 

j<k 



ife) 2 



(3.7) 



The a;-dependent density profile of counter-ions is de- 
fined as the average n(x) = ^( r — r j))- The mean- 
value calculation with the energy (|3. 71) is simplified in 
two ways. Firstly, since the y-coordinates do not mix 
with x-coordinates, they can be neglected in the consid- 
ered order. Secondly, since all (identical) particles are 
exposed to the same one-body potential of the charged 
line, a summation over particle degrees of freedom can 
be represented by just one auxiliary coordinate. We find 



n(x) = Ce x dx e" 
Jo 



E 



(x-x') 2 



27r 2 r ^ r 

j=-oo 



-o 



r 2 



(3.8) 



where C is determined by the normalization condition 
(|2.7[) . Simple algebra yields 



h(x) = c 



1 + 



(3.9) 



The contact theorem (|2.8p is fulfilled by this profile. 
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B. Two charged lines 



where C is determined by the normalization condition 



In the problem of two lines at distance d, carrying the 
same charge density ere, the electric field between the 
lines vanishes. At zero temperature T — > oo, the counter- 
ions collapse on the lines, see Fig. 2b. The Wigner crys- 
tal is thus composed of two one-dimensional arrays of 
sites with the lattice constant a = l/er, shifted with re- 
spect to one another by half-period a/2. The vertices 
of the Wigner lattice will be denoted as — (0, j/a) 
(j = 0,±1,±2, . . .) if they belong to the line at x = 

and as R^ = (d,j/a) (j — ±|,±|, . . .) if they belong 
to the line at x = d. The SC regime of large T lends 
itself to analytic progress when the inequality d <C a, or 
equivalently 



dx n(x) = 2. 



After simple algebra, we arrive at 




(3.14) 



O ^ 



da < 1 



(3.10) 



(3.15) 

This expression has the needed x — > d — x symmetry. 

To derive the pressure between the lines, we apply the 
contact theorem (|2.1[) to obtain, in the SC T — > oo limit, 



(3.16) 



- 2 2d 
P =-~-l + — 

d 9r 




is obeyed. As before, the shifts of particles from their 
Wigner positions along the y direction have no effects on 
statistical averages in the leading SC order and in the first 
1/r correction, so we shall not consider them. Let first 
one of the particles, say the one localized on the x = 
line at Rq '' = (0, 0), be shifted along the a;-direction by a 
small amount x > 0. The corresponding energy change, 
up to harmonic terms, is given by 



j38E(x) ~ — x? - — \d 2 - (d- x) 2 



(3.11) 



where we used the formula 3 l/j 2 — n 2 12. 

'3— 3>5> — 

When all particles are shifted in the x direction within 
the area limited by the two lines, R^ — > (xf\j/a) (j — 

0,±1,±2,...) andRW -> (xf\j/a) (j = ±\, ±§, . . .), 
the energy change is given by 

1 ■ 



PSE({xj}) 



4 0) ] 2 



27r 2 r ^ 



(j - W 



i 



\x (d) 



4 d) ? 



2n 2 T^ 

3<k 



1 



2vr 2 r 



E 

3 1* 



U - k f- 



(J - kf 



(3.12) 



The density of counter-ions at x (0 < x < d) is calcu- 
lated with the energy (|3.12l) : 



n{x) 



C / dx' 
Jo 



1 



2tt 2 T 



E 



(x - Z' f 



2vr 2 r 



E 



j=±l,±2,... 

(x - x'f 

J 2 



O 



r 2 



(3.13) 



Note that the small-d pressure expression (|2.14p . ob- 
tained in the weak-coupling limit, and (|3.16p obtained in 
the strong-coupling limit, coincide only in the leading 2/d 
term. The sub-leading constant terms differ. It should be 
kept in mind here that the present expansion makes sense 
provided Eq. (|3. 101) is fulfilled, which means that d -C V. 
Given that Q3.16P exhibits an attractive regime, for large 
enough V, for values of d slightly above 2, it can be con- 
cluded that Eq. (13. 16)) is able to capture the possibility of 
like-charge attraction, a phenomenon first reported some 
30 years ago [43Tj45| . It is also interesting to comment on 
the mean-field failure to capture such an effect, in spite 
of a small d expansion of the pressure, Eq. (|2.14j) . that is 
very close to its SC counterpart p,16j) . The requirement 
for the validity of expansion (|2.14l) is that d -C 1, out- 
side the regime of d > 6 where Eq. (|2.14l) , taken without 
particular precaution, would lead to a negative pressure. 
The requirement d <^ 1 stands in the strongly confined 
regime, where the entropic cost of confining the counter- 
ions is overwhelming, and leads to a dominant repulsive 
pressure P ~ 2/d. 

Taking the expansion order as indicated in formula 
([3T6| . the attractive (P < 0) and repulsive (P > 0) 
regions in the (r, d) plane are split by the solid curve in 
Fig. 3. The maximum attraction, given by dP/dd = 0, is 
obtained for d max = 3-s/T, see the dashed line. Open sym- 
bols for T = 4, 6 are the results obtained in the next sec- 
tion from an analysis of finite particle numbers. It is im- 
portant here to emphasize that the equation of state lead- 
ing to Fig. 3 is trustworthy for d -C T. As a consequence, 
the dashed line showing the locus of maximal attraction 
becomes asymptotically exact (since d max /T oc r -1 / 2 
vanishes at large coulombic couplings), but the upper 
part of that diagram in Fig 3 is only indicative, as the 
question mark indicates. It will in particular be shown 
in the next section that the re-entrance of the repulsive 
regime at fixed T, i.e. the fact there exist two distances d 
where the pressure vanishes, does not take place at T = 4 
nor at T — 6. For such couplings, the first zero of P 
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defined as 



FIG. 3. Phase boundaries following from the equation of state 
(|3.16p . The solid curve corresponds to P = 0, and discrim- 
inates the repulsive region on the left hand side from the 
attractive one. The dashed curve corresponds to the max- 
imum attraction, given by dP/dd = 0. The dashed-dotted 
line stands for the phase boundary following from the refined 
equation of state (|4.47p . where the function of coupling Am is 
given by its infinite N extrapolation (|4.50[) . The two circles 
stand for the exact location of the P = points for Y = 4 
and 6, obtained in section 4. Likewise, the two squares are 
the corresponding exact maximum attraction points at V — 4 
and 6 (see text). We stress that such diagrams are obtained 
from strong coupling in conjunction with short distance ex- 
pansions. Consequently, for the range of T values plotted, a 
qualitative picture only can be expected. The question mark 
is a reminder of the fact that d should be much smaller than 
T to allow for our expansion. 



only (with d slightly above 2) is observed. Consequently, 
while the re-entrance phenomenon is a feature backed by 
Monte Carlo simulations for three dimensional systems 
[l7| . its existence in the present two dimensional situa- 
tion is still an open question, that is difficult to address 
with analytical tools. 



1 f N 

Zn = n\ / n i d2z 3 n 

' JA j=l 3 <k 



Zk\ 



(4.1) 



where w(z, z) = exp[T^(z, z)] is the one-body Boltzmann 
factor. In Zn is the generator for the particle density in 
the following sense 



n(z, z) = w(z, z) 



8w{z, z) 



In Zn- 



(4.2) 



For the couplings T = 2j with 7 a positive integer, a 
ID Grassmann representation of the 2D partition func- 
tion (|4.1[) was derived and further developed in a se- 
ries of works [32TI341 ] . Let us introduce on a discrete 
chain of N sites j — 0, 1, . . . , N — 1 two sets of Grass- 
mann variables {^ Q \ each with 7 components 
a = 1, ... ,7. The Grassmann variables satisfy the or- 
dinary anti-commuting algebra [46j . The partition func- 
tion (|4.1[) is expressible as an integral over the Grassmann 
variables in the following way 



S(£,V0= E E '">^- 

j,fc=0 



(4.3) 



Here, VipT>£ = UfJo'^j 

the action S involves pair interactions of composite op- 
erators 



(7) 



,dd 1} and 



E 

■ ■ ■ ■ +i-> 

E 



(7) 



3l >••• J-y=0 

(j'i + ---+j7=i) 



Ai) 

' Sfc., ) 



(4.4) 



k 1 , . . . , k-y —0 

(fei+...+fe 7 =fe) 



i.e. the products of all 7 anti-commuting-field compo- 
nents, belonging to either £-set or ip-set, with the fixed 
sum of site indices. The interaction matrix has the ele- 
ments 



IV. "INTERMEDIARY" FLUID REGIME 



General formalism 



We consider N identical particles of charge — e, con- 
strained to a 2D domain A; points of the domain will 
be specified by the complex coordinates z = x + iy, z = 
x — iy. The particles interact pair-wisely through the 
2D Coulomb potential v(z, z) = — In \z — z'\ and are ex- 
posed to a one-body potential eip{z, z) due to the uniform 
charge density ere on the line walls forming the domain 
boundary dA. The partition function at coupling T is 



w jk = d 2 zw(z,z)z^z k ; j,fc = 0,l,...,7(JV-l). 

J A 

(4.5) 

The representation (|4.3[) provides Zn(j) as a function of 
interaction elements {wjk}- The particle density (|4.2p is 
given by 



■y(N-l) 

n(z,z) = w(z,z) ^2 



k)z j z k , 



where the two-correlators 



Z N 




dw 



r{V J 

Id.Z n (j). 



(4.6) 



(4.7) 
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The above Grassmann formalism is straightforwardly 
applicable to the case of the cylinder surface with the 
Coulomb potential (|1.3|) . Due to the periodicity along the 
y axis, both the one-body Boltzmann factor w and the 
particle density n are only ir-dependent. The interaction 
Boltzmann weight for two particles at the points z = 
x + \y and z' = x' + \y' is expressible as 



2 sinh 



7r(z — Z 1 ) 



W 



_ e -nV{x+x')/W 



e -2irz/W _ e -2irz' /W 



(4.8) 



For each particle of the iV-particle system, the prefactors 
from N — 1 pairwise interaction Boltzmann weights and 
the multiplication by the constant An/W 2 (which is ir- 
relevant from the point of view of the particle density) 
renormalize the one-body w{x) in the following way 



Wrcn(x) 



All 

w 



2 w(x) exp 



2ttj 
~W 



(N- 1)5 



(4.9) 



The partition function is again given by (|4.1[) . with the 
substitutions w — > w lcn and z — > exp(— 2ttz/W). Due to 
the orthogonality relation 



dy exp 



= ws jk , 



(4.10) 



the interaction matrix (|4.5j) becomes diagonal, Wjk — 
WjSjk with 

Wj = W J dxw rcn (x)exp (-^jx^j . (4.11) 
Due to the "diagonalized" form of the partition function 

f j(N-l) 

Z N { 1 )= Y[ exp&wM), (4.12) 



3=0 



only two-correlators 



(Si*i) 



d 



lnZjv(7) 



will be nonzero. The density is thus given by 

7 (W-i) 



n(x) = w T 



(x) ( H i*j>exp(-Si«) 

3=0 V J 



(4.13) 



(4.14) 



We see that the original problem reduces to finding 
the explicit dependence of Zn(h) 011 the set of weights 
{ w j}]=a 3 sa y by using the anti-commuting integral in 
(14121 . 

Zn(i) can be found trivially for 7 = 1 (T = 2) 
when the composite operators are the standard anti- 
commuting variables Sj = £j , *Sfj — ijjj . The partition 
function contains the only term 



Z N (1) = w wi • • • w N - 



(4.15) 



Consequently, 



(~ j y j } = — for all j = 0,1,..., N- 1. (4.16) 

In the case of higher integer 7's, the number of (always 
positive) terms increases quickly with N. For N = 2 
particles and arbitrary integer 7, we have 

For N — 3 particles, we have 

^3(2) = l(JoW2W4 + 2z«ow| + 2w\wi 

+4wiw 2 w 3 + 6^2, (4-18) 
Z 3 (3) = + 3 2 wou;4W5 + 3 2 wiw 2 w e 

+6 2 WiWsW5 + 15 2 W2W3W4, (4-19) 



To document the number of terms in Zn(^/) 
mention that when all w, = 1 then Zn(~/) — 



etc. 

we mention tnat wnen all Wj 

( , yN)\/[('y\) N N\\. The methods for systematic generation 
of Z]y(j), realized in practice through computer language 
Fortran, are summarized in Ref . [33| . We were able to go 
up to N = 10 particles for 7 = 2 and up to N = 9 parti- 
cles for 7 = 3. For the sake of completeness, we mention 
that the number of terms in Zn(j) is on the order of 10 8 
for 10 particles at 7 = 2. 

For 7 being an odd positive integer, the composite 
operators S and \& are products of an odd number of 
anti-commuting variables. This is why they satisfy the 
usual anti-commutation rules {Hj, S&} = {^j, ^k} = 
{SjjVE'/j} = and, in particular, we have S| = = 
0. Each exponential in (|4.12[) is then expanded as 
exp(SjWj i Sj) = 1 + 'EjWj'&j. We conclude that, for odd 
7, a given Wj can occur in a summand of Zn("j) at most 
once. In view of (|4.13[) . this property implies the inequal- 
ity 



"3(^3. 



<1 for 7 =1,3, 5,. 



(4.20) 



On the other hand, if 7 is an even positive integer, the 
composite operators are products of an even number of 
anti-commuting variables and therefore commute with 
each other: [Ej,S k ] = [*j,* fc ] = [Ej,^ k } = 0. Higher 
powers of Wj are then allowed in summands of Zjq (7) and 
the inequality (|4.20[) has no counterpart for even 7's. 



B. Cylinder: Single charged line 

We consider the periodic strip of circumference W , 
semi-inhnite in the x-direction, x € [0, 00], see Fig. lc. 
The charge density o~e at line x = is neutralized by 
N = uW particles of charge — e. The potential induced 
by the line charge is — naex, so that w{x) = exp(—T-!rax) . 
The renormalized one-body Boltzmann factor (|4.9p and 
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the interaction strengths (|4.11[) take the form 



4tt 



W *M = exp ( -—X ) , Wj 



2lT"f 



w 



1 



The particle density (|4.14l) reads 



J + (7/2)' 

(4.21) 



n(x) 



An 
W 2 



J'=0 



47T 

W 7 



+ i) 



(4.22) 

For finite N (or, equivalently, finite W), the particle den- 
sity exhibits at asymptotically large distances an expo- 
nential decay to zero, lim^-^oo n(x) oc exp(— 2ttjx/W). 

Our aim is to continualize the formula (|4. 22[) in the 
thermodynamic limit N,W — > oo, at the fixed ratio 
N/W = cr; this makes the semi-infinite cylinder sur- 
face equivalent to the system of the charged straight 
line in contact with half-space occupied by counter-ions. 
For a given finite N, we define a set of discrete values 
& = 7^-<Sj%) with j = 0,1,..., 7 (JV - 1). As the 
continuous variable, we choose t = j/[y(N — 1)], tak- 
ing values in the interval [0,1]. In the continuum limit 

TV — s- oo, the set of discrete values fy^ tends to a con- 
tinuous positive function 



= Jim /W, 



JV-i-oo 



t = i/[7(^-i)], 

i = 0,l,..., 7 (iV-l). 

(4.23) 

The continualization of (|4.22l) results in 

h(x)=2 f d«/ (7) (i)exp(-2te). (4.24) 
Jo 

We see that the original problem reduces to the one of 
finding the function f^\t). The electroneutrality con- 
dition (|2.7I) and the contact theorem (|2.8|l hold provided 
that the function f^'(t) is constrained by 



dt/W(t) = l, 



f dttf™{t) = ~, (4.25) 



respectively. 

Let us first perform a brief analysis of the general den- 
sity formula (|4.24l) . without knowing explicitly f^\t). 

For odd 7, the inequality (|4.20|) implies that /J^ < 7 for 

all j, and so < 7 in the whole interval i g [0,1]. 

Consequently, 



n(5) < 



7 
2i 2 



[l - (1 + 2£)c~ 2i ] for odd 7. (4.26) 



In particular, h(x) < j/(2x 2 ) at large x. This already 
provides a non-trivial bound. Note that this relation 
teaches us that the present analysis, valid for integer val- 
ues of 7, cannot be continualized to small couplings to 
encompass the mean field limit 7 — > 0, since at mean-field 
level, one has h ~ x~ 2 . It is furthermore clear that the 
asymptotic decay of the particle density is determined by 



the behavior of f^(t) in the limit t — > 0. Let us assume 
that /W(t) has a power-law behavior f^\t) ^ t ^a ct u , 
where the positiveness and the boundedness of /^ 7 ' (t) are 
ensured by c > and v > 0, respectively. The special 
case of z/ = corresponds to the situation when f^\t) 
approaches a positive number as t — > 0. Inserting our 
assumption into (|4.24[) . we get 



fi(x) ~ (1 + v)\ 



1 



(4.27) 



This means that if > for larger couplings, the asymp- 
totic decay of the particle density is faster than the weak- 
coupling prediction (|2.6|) . 

For 7 = 1, we have the exact result /j 1 ^ = 1 for all j = 

0. 1. . . . , N — 1 and particle numbers N, so that /W (t) 

1 . This function fulfils the normalization relations (|4.25p . 
The formula (|4.24l) for the density profile becomes 



1 



K*)=2p[l-(1 



2£)e 



-2x1 



At large i, 



x —> 00. 



(4.28) 



(4.29) 



This behavior resembles, up to the renormalization factor 
1/2, the weak-coupling decay (|2.6[) . The independence of 
the counter-ion density on the line charge density ere is 
present also at the considered finite temperature. 

For 7 = 2 and 7 = 3 (i.e. r = 4 and 6), we performed 
exact calculations up to N = 10 and N = 9 particles, 
respectively. The discrete representations {fj%}^=o ^ 
and {/j ^r}jlio _1 ^ °f tne corresponding functions f^ 2 \t) 
and f( 3 \t) are presented in Fig. 4. It is seen that data 
converge rather quickly when increasing particle num- 
bers. The continuous functions f^ 2 \t) and f^ 3 '(t) are 
well approximated by the corresponding discrete plots 
obtained for N = 10 and N — 9 particles (dashed lines). 
Although the plots look at first sight to be symmetric 
with respect to t = 1/2, they are not. 

In view of the above discussion about the large- 
distance decay of the particle density (|4.27p , the behavior 
of the functions f^ 2 \t) and f^(t) in the limit t -> is 
of primary importance. The first relevant question is 
whether / (2) (0) and / (3) (0) are positive or equal to 0. 

Since t = for j — 0, we have f^(0) = limjv-!.oo /q jv 



Given that the "natu- 



and /< 3 )(0) = lim^oo / > 
ral" variable in the finite-iV analysis is x = 1/(N — 1), 
the dependence of the sequences {/n 2 ^-} and {fo%} on 
1/(N — 1) is pictured in Fig. 5. For 7 = 2, the sequence 
is well fitted by the linear form /q ' n ~ 0.195 + 0.85a; 
with x = 1/(N - 1) (dashed line), i.e. / (2) (0) ~ 0.195 
is apparently positive and v = in (|4.27[) . We note that 
this value is fully corroborated by the analysis of the be- 
haviour of f[ 2 \ that is well fitted by f[ 2) N ~ 0.195 + 1.70a; 
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FIG. 4. Discrete representations of the functions f^(t) (V = 
2 7 = 4) and f m (t) (T = 6) for increasing number of particles 
AT. 



(dashed line in Fig. 5). We conclude that at 7 = 2, 
the large-distance behavior of the density exhibits the 
mean-field behavior 1 jx 1 , with a renormalized prefactor 
f( 2 \0)/2 ~ 0.1, significantly smaller than its mean-field 
value of 1, or the 1/2 value which holds at T = 2. This 
illustrates enhanced screening at larger Coulombic cou- 
pling r. 

For 7 = 3, the sequence is well fitted by the quadratic 
form f^ 3) N - 0.2x + 0.32a; 2 (dashed line), i.e. / (3) (0) - 
and v > in the large distance behavior (|4.27l) . The 
sequence {/q%} does not provide any information about 
the value of v. The index v can instead be deduced from 
the subsequent sequence {f[ 3 jf} with the coordinate t = 
l/[3(iV— 1)] which, as N increases, mimics the plot of 
the function /( 3 ) (t) for small t. In the limit N — > 00 (t — 

0), the sequence {/-{ ^} must converge to the previously 
obtained / (3) (0) = 0. The sequence is well fitted by a 



power-law /£ 



(3) 



c u with v ~ 1.45, see the solid line 



in Fig. 5, indicating that v ~ 1.45. The stability of the 
fit is documented in the inset, that shows the sequence 

(3") 

{f{ jy} on log-log plot; the dashed line has slope 1.45 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 H 
0.1 





(2) 



0.1 0.2 0.5 

(N-l) 



f 



(2) 




(N-l) 



FIG. 5. The jV-dependence of discrete sequences /q 7 j^, f^jtf for 
7 = 2, 3; their relation to the small-t behavior of the functions 
f^ 2 \t) and f^ 3 \t) is explained in the text. The symbols are 
for the exact results at finite N. For 7 = T/2 — 2, both f[ 2 ^ 
and /q extrapolate to the same limiting value (0.195) for 
N — ► 00 (in other words, v = in Eq. (|4.27|) . On the other 
hand, for 7 = 3, a similar analysis provides a vanishing limit, 
i.e. v 7^ 0. In the main graph, the dashed lines are linear 

fits while the continuous curve corresponds to a power-law 

(3) 

with exponent 1.45. The inset shows f{ on a log-log scale, 
to evidence power law behaviour. The dashed line has slope 
1.45 while the dotted line, given as a guide to the eye, has 
slope 1.6. 



and the dotted line has slope 1.6. A similar value of 

(3) 

v is obtained from the sequence \f\ ^}, and also from 
a complementary direct plot of the data of Fig. 3 for 
r = 6 on a log-log scale (not shown). We therefore see 
that at 7 = 3 the density behaves like ~ l/i a at large 
distances with an exponent a = 2 + v close to 3.5, in 
contrast to the mean-field prediction. For this coupling, 
the asymptotic density decay depends on the line charge 
ere. Indeed, returning to original variables, one has in 
general 

/ \ — v —2 — 1/ 

n{x) cx a x 



C. Cylinder: Two charged lines 

We now consider the periodic strip of circumference 
W, finite in the x-direction, x € [0, d], see Fig. Id. The 
equivalent charge densities ae at lines x = and x — d 
are neutralized by N = 2aW particles of charge — e. The 
electric field generated by the line charges vanishes, so 
that w(x) = 1 and 



47T 



w v 



l)x 



(4.30) 
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With respect to the equality d/W = d/(irjN), the inter- 
action strengths (|4.11[) take the form 



1 



w 



' ■ -( V-l) I 



\ _ g -yN | 



J 



-*("-!)]} (4.31) 



for j ^7(JV-l)/2 and = 4d/(jN) for j = 7 (A/-l)/2. 
The density profile is given by 



4tt 



7 (AT-1) 



E (Hi^e-^-i^K (4.32) 

The reflection symmetry n(x) = n(d — x) implies the 
relation 

w j( E j^j) = ^ 7 (W-l)-j< S 7(JV-l)-J*7(iV-l)-i) ( 4 - 33 ) 

which is valid for all j = 0, 1, . . . , j(N — 1). 

The continualization of the formula (I4.32[) . in the limit 
N, W — > oo while keeping the ratio N/W — 2a, proceeds 
along the above lines. As the continuous variable, we 
choose t = 2\j-f(N-l)/2]/[y(N- 1)], taking values in 
the interval [—1,1]. In the continuum limit, the discrete 

values fjj^{d) = jWj(Ej^Jj) tend to a positive bounded 

function f^'(t,d). The continualization of (|4.32l) results 



dii/ (7) (i, d)- 



h{x) — 

The symmetry n(x) = n(d — x) implies 



(4.34) 



(4.35) 



which enables us to rewrite (|4.34[) into a more convenient 
form 



n{x) =2 dttf { ~>\t,d)- 



cosh 



2/ 



M) 



sinh(ic?) 



(4.36) 



The electroneutrality condition (|3.14[) leads to the con- 
straint 



dt/^(M) = i. 



(4.37) 



According to the contact theorem (|2.ip . the (renormal- 
ized) pressure is given by 



i.e. 



P = 2 I dtt 
Jo 



P = n(0) - 1, 



f {n \t,d) coth(id) - 1 



(4.38) 



(4.39) 



For r = 2, we have the exact result /ijy-(d) = 1 for 
all j = 0, 1, . . . , N — 1 and particle numbers N, which 



implies f^(t, d) = 1 in the whole interval t € [0, 1]. The 
pressure 



P = 2 



exp(— td) 



dtt- 
II smh(td) 



(4.40) 



is positive for every distance d, so there is always the 
repulsion between two equivalently charged lines at T — 
2. For small distances d, we have 



P = i 
d 



2d m i . 



(4.41) 



This expansion resembles the SC one (|3T6]) , up to the 
renormalization factor 2 ahead of the d term. For large 
distances between the lines, the formula (14.40[) yields 



d^ 



(4.42) 



This result coincides, up to a renormalized pref actor, 
with the PB prediction (|2.15j) . The asymptotic decay 
of the pressure is universal in the sense that it does not 
depend on the magnitude of the charge density ae on the 
lines. 

It is instructive to compare the exact solution (|4.40[l . 
obtained in the thermodynamic N — > oo limit, with the 
finite- TV results. For finite N, equations P~52|) and 
imply 



i(N-l) 



(4.43) 



For 7 = 1, we have (Ej\Pj) = 1/wj for all j = 
0, 1, . . . , N — 1. To derive the asymptotic d — > oo be- 
havior of P, we note from (|4.3ip that for 7 = 1 



Jim — 

d— >oo 



j-\(N-l) if j-l(N-l) >0, 
if j- |(JV — 1) <0. 



(4.44) 

Consequently, for even particle numbers iV, we have 



Jim P 



iV 2 



N-l 

■ E 

j=iV/2 



1 



i--(^v-i) 



1 = 0, (4.45) 



while for odd TV: 



8 N ~ 1 

j = (W+l)/2 



1 



1 

iV2- 



(4.46) 

The above results can be understood intuitively by the 
discreteness of particles. If the number of particles is 
even, i.e. N = 2N* , at asymptotically large distance 
each of the charged lines attracts just N* particles. The 
whole system thus consists of two neutral subsystems 
which do not interact with one another. On the other 
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hand, if the number of particles is odd, N = 2N* + 1, 
one of the particles is shared by both lines, so that the 
two double-layers can never strictly decouple. This "mis- 
fit" particle is responsible for the asymptotic attraction 
— 1/N 2 between the charged lines in (|4.46[) . The asymp- 
totic attraction disappears in the thermodynamic limit 
N — y oo. It is interesting that exactly same asymptotic 
relations (|4.45l) and (|4.46l) can be found for higher values 
of 7, as was verified on finite- N calculations for 7 = 2,3. 
A detailed discussion about this interesting finite- N phe- 
nomenon will be published elsewhere (47j . 




FIG. 6. Pressure curves for F = 2. In the main graph, the 
pressure is shifted by 1/N 2 for odd values of particle numbers 
N. For the asymptotic N — 5- oo result shown by the thick 
line, see formula (|4.40l) . 

The r = 2 results for the pressure dependence on the 
distance between two lines in the case of finite N are 
compared with the asymptotic N — > oo result (|4.40[) in 
Fig. 6. In the main graph, the pressure is shifted by 
1/N 2 for odd values of N (so that it vanishes at large 
distances) while curves with even N are not shifted. The 
decay of the pressure is always monotonous. We see that 
curves with even TV lie below and the shifted ones with 
odd N lie above the asymptotic N — > oo line; the odd 
and even curves systematically "sandwich" the asymp- 
totic one. The inset, where only even TV are considered, 
shows a quick convergence of the pressure plots for finite 
particle numbers to the asymptotic one. To our surprise, 
relatively small numbers of particles are sufficient to have 
realistic estimates of the pressure in the thermodynamic 
limit, at least up to "reasonable" distances d ~ 5. 

The results for the pressure dependence on the distance 
in the case of the smallest particle number N = 2 and F = 
2, 4, . . . , 12 are depicted in Fig. 7. For r = 2, we have 
the already discussed monotonous decay of the pressure 
to 0. For T > 4, the pressure becomes negative at a In- 
dependent distance and remains negative up to d — > oo, 
i.e. there is no further intersection of the curve with 
the P = axis. We learn that the attraction between 
equivalently charged lines is not associated only with the 
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FIG. 7. Rescaled pressure versus distance between two 
charged lines for N = 2 particles and the couplings F — 
2, 4,..., 12. 



thermodynamic limit, but manifests itself (at sufficiently 
large T) even for N = 2 particles. 

Pressure curves for T = 4 and T = 6 in Fig. 8 are 
presented separately for even N — 2, 4, 8 even (the main 
graph) and odd N = 3,5,7 (the inset). We see that 
by increasing N the difference between plots becomes 
very small and so they are probably very close to the 
asymptotic line, at least for d < 5. As in the N — 2 case, 
the pressure crosses the P = line at some distance. This 
distance, evaluated at N — 8, is indicated for r = 4, 6 
by open circles in Fig. 3; we see a good agreement with 
the SC phase diagram (solid line). On the other hand, 
the distance of the maximum attraction (open squares) 
is relatively far away from the corresponding dashed line; 
this might be caused by an extended plateau around the 
minimum point. As before, after crossing the P = 
line, the curves remain in the attraction region up to 
d — > oo. This fact sheds doubts on the existence of the 
upper branch of the phase diagram (for N — > oo) in Fig. 
3. We cannot answer this question because the finite- N 
calculations do not reflect adequately the pressure in the 
thermodynamic limit just in the region of large d. 

We end up this section by an analysis of the small- 
distance expansion of the pressure; see formula (|2.14l) for 
the weak-coupling regime, (|3.16[) for the SC regime and 
(I4.4ip for T = 2. The leading small-distance term is 
always 2/d as a consequence of the spatial homogeneity 
of the particle density h(x) ~ 2/d in the limit d — > 0. 
The next (constant) term is equal to —1/3 in the weak- 
coupling limit and to —1 in the SC limit and for F = 2. 
The value —1 was detected also in finite- A calculations 
for T = 4, 6 -as Fig. 9 shows- and is expected to persist 
up to r — ¥ oo. From the point of view of the contact 
theorem (|4.38p . this means that in the limit r — > the 
contact density behaves like n(0) ~ 2/d + 2/3 + 0(d), 
while for V = 2,4,6,... it behaves like n(0) ~ 2/d + 
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FIG. 8. Pressure curves for r = 4 and T = 6. Results are 
presented separately for even TV = 2, 4, 8 (the main graph) 
and odd N = 3, 5, 7 (the inset). 



It is quite remarkable that this simple form accounts for 
the exactly know results at T = 2 and r — > oo, to- 
gether with the numerically accurate data in the ther- 
modynamic limit obtained at T — 4 and T = 6. It is 
however inapplicable to the mean- held limit r — > 0. 
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0(d). It is not clear at which T < 2 the fundamental 
change in the short-distance behavior of n(0) starts; in 
the subsequent analysis, we shall restrict ourselves to the 
region T > 2, where 

P= ~- 1 + A N {T)d + 0{d 2 ). (4.47) 
d 

We have at our disposal the exact information about the 
r — > oo limit of the prefactor: 



9r 

lim — Aoo(r) = 1. 
r— >oo z 



(4.48) 



The simplest Pade approximant corresponding to this 
asymptotic formula is 



4»(r) 



2 r 



9r r + b 



(4.49) 



where the free coefficients a, b should be fixed by the 
known (exact or approximate) values of Aoo(r) at some 
T-points. We know the exact result at T = 2: = 
2) = 2/9. We were able to evaluate accurately A^T) 
for r = 4, 6; the analysis for T = 6 is presented in 
Fig. 9. Considering P — 2/d + 1 as the function of 
d in the main graph, An is nothing but the slope y/x 
of the N — 6, 7, 8 data curves, taken in the small- a; 
limit. In the inset, the slopes An (N — 6,7,8) are fit- 
ted against l/N with the function y = 0.0454 - 0.013a;, 
hence the value AoofT = 6) ~ 0.045 for large N. The 
same scenario applies to T — 4, where the limiting slope 
is A^r = 4) ~ 0.075. The three results at T = 2,4,6 
are perfectly matched by the two-parameter Pade ap- 
proximant (I4.49P if we choose a = 1/4 and b = —7/8, 
i.e. 



2 r + i/4 

9T T-7/8' 



r > 2. 



(4.50) 



FIG. 9. The procedure of fitting the coefficient ^^(r = 6), 
defined by the expansion (|4.47|l . A jv corresponds to the slopes 
of the JV = 6, 7, 8 data sets in the main graph. In the inset, An 
are fitted linearly against l/N to obtain Aoc(T = 6) ~ 0.045. 
Note the vertical scale in the inset that indicates that the 
thermodynamic limit is already closely approached by systems 
with N = 6 particles. Very similar results were found at 
r = 4, with different numerical constants. 



The result (|4.50|) . when substituted into the expan- 
sion (I4.47P , can be used to improve the phase diagram in 
Fig. 3, that follows from the SC equation of state (|3. 161) . 
The problem of our original SC phase diagram of Fig. 
3 is that it includes T = 2 in the attractive regime (for 
some distances d around 4), while we have shown that 
for N —¥ (X), there is no attraction between lines at the 
exactly solvable T = 2 case (see e.g. Figure 6). It can be 
seen in Fig. 3 that making use of the improved equation 
of state (|4.47|) together with (|4.50[) . leads to a shift of 
the attraction border towards higher couplings (see the 
dashed-dotted line). More precisely, the critical T below 
which no attraction is possible is changed from a value 
1.77 with to 2.81 with P~5D]) . Since our exact re- 

sults indicate the possibility of attraction at V — 4 and 
6, but none at T = 2, the second estimation seems more 
reliable. 



V. CONCLUSION 

In this paper, we have studied 2D models of counter- 
ions at and between charged lines. Since only one type of 
counter-ions was considered (no salt), these ions can be 
considered as point-like without any ensuing pathology. 
2D models, while maintaining the essence of 3D Coulomb 
models, are simpler to handle analytically as well as nu- 
merically. In the one-line geometry, we focused on the 
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large-distance decay of the particle density. In the two- 
line geometry, the small distance behavior of the pressure 
was the center of interest (although of particular signif- 
icance, the large distance behaviour is more difficult to 
obtain, and could only be addressed in particular cases). 
The possibility of an attraction (negative pressure) for 
some distances between equivalently charged lines, medi- 
ated by counterions, was investigated. 

The weak-coupling limit has basically the same mean- 
field Poisson-Boltzmann (PB) form in any dimension. 
For one-line geometry, the counter-ion density falls at 
asymptotically large distances like the inverse-power-law 
1/x 2 which does not depend on the magnitude of the line 
charge ae. The same phenomenon is observed in the two- 
line geometry: The asymptotic decay of the pressure does 
not depend on ae. This phenomenon occurs in both 2D 
and 3D. The pressure is always positive: the attraction 
phenomenon is absent for weak couplings. 

The strong coupling (SC) anal ysis presented here is 
the 2D adaptation of the method 27J based on the har- 
monic expansion of the interaction energy around the 
ground-state Wigner crystal formed by counter-ions. The 
method is applicable to small distances. The small- 
distance expansion of the pressure between two lines pro- 
vides a phase diagram -in Fig. 3- which includes the 
attraction region. Our strong-coupling expansion dif- 
fers from those that have previously been proposed: it 
allows to compute the corrections to the leading order 
term, which is not the case of the original method of 
Netz and collaborators 0, [lfl [l|| , that nevertheless suc- 
cessfully predicts the leading term. Our approach also 
differs from that of Refs. [l8|, [llj], where the excitations 
considered (counter-ions displacements restricted to the 
charged interface) are not those that turn relevant at 
short distances. As a consequence, the predictions of 
Refs. [HI, [lj| do not cover the short-range phenomena 
that we have studied here under strong coupling. 

The intermediary (i.e. between weak- and strong- 
coupling regimes) values of T — 2,4,6 are studied by 
using a previously developed Grassmann variables for- 
malism [32-34]. As a constraining domain for counter- 
ions, we choose the surface of a cylinder; this enables us 
to mimic infinite systems by finite ones containing finite 
numbers of particles N. We systematically observed that 
a system with as little as N = 5 to 10 particles may be 
considered as "large" , in that it is already close to the 
thermodynamic limit. 



The case r = 2 is solved exactly, for finite as well as 
infinite TV. It shares many features with the PB theory: 
The particle density decay in 1/x 2 does not depend on 
ae, the pressure is always repulsive and its large-distance 
asymptotic is independent of ae. 

The couplings T = 4, 6 are investigated for finite par- 
ticle numbers N. In the one-line problem, the particle 
density decay is still of the PB type 1/x 2 for T = 4. For 
r = 6, there are strong indications that the counter- ion 
density behaves like 1 / x a for large x where the exponent 
a is close to 1.5 (we found 1.4 < a < 1.6). This sig- 
nals the breaking down of the large-distance PB theory 
and the dependence of the asymptotic density on ae. It 
has been argued on the contrary that a strongly coupled 
double-layer behaves, at large distances where counter- 
ions correlations should be less important, as predicted 
by a suitably renormalized mean-field approach [3^, H|| . 
The data reported here provide evidence, for two dimen- 
sional systems, that this is not always the case. The 
corresponding question for 3D systems that are the main 
objects of interest in Refs [HI, [36| (i.e. with 1/r inter- 
actions instead of logr) remains open. As concerns the 
pressure between equivalently charged lines, it is highly 
non-trivial even for N — 2 particles (see Fig. 7). The 
monotonous decay at T = 2 changes for T > 4 and turns 
into a profile with an attractive regime, starting from 
a certain distance. Increasing N, the results converge 
quickly for intermediate distances of interest. We see 
that the attraction phenomenon is not restricted to the 
thermodynamic limit, but takes place for small particle 
numbers and relatively small values of T (4 and 6). The 
finite- N analysis of the equation of state enabled us to 
improve the phase diagram evaluated in the SC regime. 

Many questions are still open. Among interesting per- 
spectives is the question the effective interaction between 
arbitrarily shaped objects. Another relevant problem lies 
in the generalization of the present ideas and methods to 
system where not only one type of micro-ions is present, 
such as electrolytes. 
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